Hamiltonian structure of the teleparallel formulation of GR 



M. Blagojevic* and I. A. Nikolic 

Institute of Physics, P.O.Box 57, 11001 Belgrade, Yugoslavia 
(February 01, 2000) 



Abstract 

O 

■ We apply Dirac's Hamiltonian approach to study the canonical structure of 

■ the teleparallel form of general relativity without matter fields. It is shown, 
r*i . without any gauge fixing, that the Hamiltonian has the generalized Dirac- 

ADM form, and constraints satisfy all the consistency requirements. The set 
of constraints involves some extra first class constraints, which are used to find 
additional gauge symmetries and clarify the gauge structure of the theory. 



(N 



I. INTRODUCTION 



> 
(N 
(N 
O 

Among various attempts to overcome the problems of quantization and the existence 
of singular solutions in Einstein's general relativity (GR), gauge theories of gravity are 
especially attractive, as they are based on the concept of gauge symmetry which has been 
very successful in the foundation of other fundamental interactions. The importance of 
the Poincare symmetry in particle physics leads one to consider the Poincare gauge theory 
(PGT) as a natural framework for description of the gravitational phenomena [0-0] (for 
^ ; more general attempts, see 0). 

Basic gravitational variables in PGT are tetrad field b ^ and Lorentz connection A 1 -^, 
which are associated to the translation and Lorentz subgroups of the Poincare group, respec- 
tively. These gauge fields are coupled to the energy-momentum and spin of matter fields, 
and their field strengths are geometrically identified with the torsion and the curvature: 
T\ v = d^V+^V^-M' Ri \v = d^v+AK^v-iiiv). The spacetime of PGT turns 
out to be Riemann-Cartan space C/4, equipped with metric and linear, metric compatible 
connection. Dynamical content of PGT is determined by the Lagrangian C = b(Cc + Cm), 
where the gravitational part Cq is usually assumed to be at most quadratic in field strengths, 
and Cm describes minimally coupled matter fields. 

General geometric arena of PGT, the Riemann-Cartan space C/4, may be a priori re- 
stricted by imposing certain conditions on the curvature and the torsion. Thus, Einstein's 
GR is defined in Riemann space V4, which is obtained from C/4 by the requirement of van- 
ishing torsion. Another interesting limit of PGT is teleparallel or Weitzenbock geometry T4, 
defined by the requirement of vanishing curvature: 
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izV(^) = °- (i-i) 

The vanishing of curvature means that parallel transport is path independent, hence we 
have an absolute parallelism. The teleparallel geometry is, in a sense, complementary to 
Riemannian: curvature vanishes, and torsion remains to characterize the parallel transport. 

Of particular importance for the physical interpretation of the teleparallel geometry is 
the fact that there is a one-parameter family of teleparallel Lagrangians which is empirically 
equivalent to GR |5],|7|||- For the parameter value B = 1/2 the Lagrangian of the theory 
coincides, modulo a four-divergence, with the Einstein-Hilbert Lagrangian, and defines the 
teleparallel form of GR, GR11. 

The teleparallel description of gravity has been one of the most promising alternatives 
to GR. However, analyzing this theory Kopczyhsky |§ found a hidden gauge symmetry, and 
concluded that the torsion evolution is not completely determined by the field equations. 
Assuming, then, that the torsion should be a measurable physical quantity, he argued that 
this theory is internally inconsistent. Hayashi and Shirafuji |TU| tried to avoid this problem 
by interpreting certain different torsion configurations as physically equivalent, i.e. related 
to each other by a gauge transformation, but the consistency of this idea in the interact- 



ing theory seems to be questionable for non-scalar matter P,ll|. Various modifications of 



the one-parameter teleparallel theory are proposed in order to avoid the above problems 



T2]JT3|] . Trying to re-examine the gauge structure of the one-parameter teleparallel geom- 



etry Nester [14] improved the arguments of Kopczyhsky ||; the predictability problem was 



stated more precisely and bound to certain special solutions. 

Hecht et al. Jl5]| traced the appearance of non-physical modes of torsion back to some 



symmetries which are necessarily present in the (3+1) decomposition of spacetime. Using 
certain geometric arguments they concluded that some components of the tetrad velocity 
are not suited to represent dynamical degrees of freedom. In other words, these velocities 
must not appear in the evolution equations, hence they should appear at most linear in the 
Lagrangian. The choice of parameters in the teleparallel Lagrangian that ensures this to 
happen is just the one corresponding to GRy. 

The teleparallel geometry possesses many salient features. Thus, Nester |16| succeeded 



in formulating a pure tensorial proof of the positivity of total energy for Einstein's theory in 
terms of the teleparallel geometry. He found that special gauge features of GRy, which are 
usually considered to be problematic, are quite beneficial for this purpose. Mielke yT7| used 
the teleparallel geometry of GRy to give a transparent description of Ashtekar's complex 
variables, while Andrade et al. |Tj| formulated a five-dimensional teleparallel equivalent of 



Kaluza-Klein theory. There are also some attempts to understand the role of torsion at the 
quantum level [19| . 



The purpose of this paper is to investigate the canonical structure of GRy using Dirac's 
Hamiltonian approach |20 , as this is, in our opinion, the best way to clarify both the nature 



of somewhat mysterious extra gauge symmetries, and the question of consistency of GRy. 
We shall find that a specific choice of coupling constants in the teleparallel Lagrangian leads 
to the appearance of some additional first class constraints, and, consequently, to extra gauge 
symmetries, which clarify the meaning of non-dynamical torsion components and give us a 
complete picture of the gauge structure of GRy. 

We remark here that Maluf [[T[|] tried to analyze some aspects of the Hamiltonian struc- 
ture of GRy. However, his approach is based on some unnecessary gauge fixing conditions, 
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adopted at the level of Lagrangian in order to simplify the calculations, so that many specific 
gauge features of the theory remained effectively hidden. 

The layout of the paper is as follows. After recalling some basic elements of the La- 
grangian teleparallel formulation of GR in Sec. II, we work out all the primary constraints 
and construct the corresponding Hamiltonian density in Sec. III. It is shown that a specific 
choice of parameters in the Lagrangian leads to additional primary constraints. Then, we 
study the consistency conditions in Sec. IV, and derive the algebra of constraints in Sec. V. 
These results are used in Sec. VI to construct extra gauge generators and clarify the nature 
of the related gauge symmetries. Section VII is devoted to concluding remarks, while some 
technical details are presented in the Appendix. 

Our conventions are the same as in Ref . f22[ : the Latin indices refer to the local Lorentz 
frame, whereas the Greek indices refer to the coordinate frame; the first letters of both alpha- 
bets (a, b, c, a, j3, 7, ...) run over 1,2,3, and the middle alphabet letters (i, j, k, ...; /1, v, A, ...) 
run over 0,1,2,3; 77^ = diag (+, — , — , — ), e 0123 = +1 and 5 = 5(x — x'). 



II. THE TELEPARALLEL FORMULATION OF GR 

Lagrangian. Gravitational field in the framework of the teleparallel geometry in PGT 
is described by the tetrad b k M and Lorentz connection A 13 '^, subject to the condition of 
vanishing curvature (|1 . 1|) . We shall consider here the gravitational dynamics determined by 
a class of Lagrangians quadratic in the torsion 0,0)11 

C = b(C T + \rR l \u + C M ) , 

C T = a{AT ljk T^ k + BT l]k T^ k + CT k T k ) = (5 l]k {T)T^ k , (2.1) 



where Xij^ u are Lagrange multipliers introduced to ensure the teleparallelism condition ( |1 . 1|) 
in the variational formalism, a = 1 /2k (k = Einstein's gravitational constant), T k = T m mk , 
and Cm is the Lagrangian of matter fields. The explicit form of /3ij k is 

Pijk = a(ATij k + BTy ik } + Cr)i[jT k ) . 

The parameters A, B, C in the Lagrangian should be determined on physical grounds, 
so as to obtain a consistent theory which could describe all the known gravitational ex- 
periments. If we require that the theory ( |2.1p gives the same results as GR in the linear, 
weak-field approximation, we can restrict our considerations to the one-parameter family 
of Lagrangians, defined by the conditions [§,0|§ 

t) 2A + B + C = 0, C = -l. 

This family represents a viable gravitational theory for macroscopic, spinless matter, em- 
pirically indistinguishable from GR. Von der Heyde p3| and Hehl pj have given certain 



theoretical arguments in favor of the choice B = 0. There is, however, another, particularly 
interesting choice determined by the requirement 

ii) 2A - B = 0, 
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It leads effectively to the Einstein-Hilbert Lagrangian £gr = —abR(A), defined in Riemann 
spacetime V4 with Levi-Civita connection A = A, via the geometric identity (|A~3T): 

bR(A) = bR(A) + b(\T l]k T> k + \T ijk T* ik - T k T k ) - 2d v {bT v ) . 

Indeed, in Weitzenbock spacetime the above identity in conjunction with the condition ( |1 . 1| ) 
implies that the torsion Lagrangian in (|2.1|) is equivalent to the Einstein-Hilbert Lagrangian, 
up to a four-divergence, provided that 

A = \, B = l, C = -l, (2.2) 

which coincides with the conditions i) and it) given above. 

The theory defined by equations (|2.1|) and (|2.2|) is called the teleparallel formulation of 
GR (GRy). Note that the equivalence with GR holds certainly for scalar matter, while the 
gravitational couplings to spinning matter fields in T 4 and V4 are in general different. 



Field equations. By varying the Lagrangian ( j2.1| ) with respect to 6*^, A %3 ^ and A 



ij 



we obtain the gravitational field equations: 

- 4V M W) - Ab/3 nm "T nmi + h»b£ T = t\ , (2.3a) 

-86%]" + 4V M (&V") = » ( 2 - 3b ) 

= > ( 2 - 3c ) 

where r 1 ^ and c"^- are the energy-momentum and spin currents of matter fields, respectively. 
The first field equation can be rewritten as 

-4V M (&A"*) + 2bf3 imn T kmn - 4b(3 nmk T nm + 5 k bC T = r k t . 



Then, combined with the identity (|Aq), it takes the form of Einstein's field equation: 

R ik {A) - \rf k R(A) = r kl /2ab . (2.4a) 

Here, on the left hand side we have Einstein's tensor of GR, which is a symmetric tensor. 
Therefore, the dynamical energy-momentum tensor must be also symmetric, T tk = T kl . 
Using the identity (XT) the second field equation can be written in the form 

V M (2a#^ + 46A,/") = <t%- , (2.4b) 

where H"? = bQifh^ — hj u h^). The integrability condition for this equation is identically 
satisfied, because the covariant divergence of the left hand side vanishes on account of 
R 13 iiv = 0, whereas 

V u a u ij —nj- Tji = . 

The first equality in this relation is a consequence of the covariant conservation of angular 
momentum for matter fields (which holds when matter field equation is satisfied), and the 
vanishing of tui follows from the first field equation. 

Simple counting shows that the number of independent field equations (|2.4b ) is 24 — 
6 = 18. The multipliers Xij af3 remain arbitrary functions of time, as will be shown in the 
forthcoming Hamiltonian analysis. For any specific choice of Xij al3 (gauge fixing), equations 
( [2.4b ) can be used to determine (at least locally) the remaining 18 multipliers \ij 0a - 

In what follows we shall investigate the Hamiltonian structure and gauge properties of 
GRy without matter fields (<fij = = 0). We expect that the results obtained here will 
be also useful for the analysis of interacting GR11. 
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III. PRIMARY CONSTRAINTS AND HAMILTONIAN 



1. The basic Lagrangian dynamical variables of our theory are Ay MV ), and the 

corresponding momenta are denoted by (n^, ir^, ^ ^u) ■ Due to the fact that the torsion and 
the curvature do not involve velocities b k and A^q, one immediately obtains the following 
set of the so-called sure primary constraints: 

fc ° = 7r fe °^O, ^ E7r y »0. (3.1) 

Similarly, the absence of the time derivative of Aj,^ implies 



)1U 



ttV«0- (3-2) 



The next set of constraints follows from the linearity of the curvature in A l \: 

^/^V-^Ai^wO. (3.3) 

Before we continue, it is convenient introduce the so-called (3+1) decomposition of 
spacetime ||22||. If n is the unit normal to the hypersurface S : x° = const., with n k = 



h>k / V*? 00 ? the four vectors (n, e a ) define the ADM basis of tangent vector fields. Introducing 
the projectors on n and S , (P±)l = n> l rik, (P\\)l — $1 ~ n * n k, any tangent vector V can 
be decomposed in terms of its parallel and orthogonal components: Vk — + nkV±, where 
V± = b k Vk, Vk = Vk — n k V±, and n k V k = 0. Using an analogous decomposition of the torsion 
and the curvature in last two indices, 

T l m k — P % fhk + 2ri[ m T J j_^ , R l3 mk = R^fhk + ^ n [mR %:> l_k] > 

we find that the parallel components T 1 ^ and R^kl are independent of velocities. The 
replacement in the gravitational Lagrangian yields C = £(T%, R^^T 1 ^, R l i ±f,n k ). 

The decomposition of eo in the ADM basis yields eo = Nn + N a e a , where N = ntb k o 
and N a = hi a b k are lapse and shift functions, respectively. We note also that b satisfies 
the factorization property b = NJ, where J does not depend on b k Q. 



Now, we turn our attention to the remaining momenta ir" ||22|| . The relations defining 
7Tj a can be written in the form 

it* = J^L- = AJ^iT) , 
9T l ±k 

where Ti k = iti a b k a are conveniently defined "parallel" gravitational momenta. Using now 
the fact that (3 is a linear function of T we can make the expansion (3(T) = /3(0) + (3(1), 
where (3(0) does not depend on "velocities" T'jj and (3(1) is linear in them, and rewrite the 
above equation in the form 

Pfk = n ik -/J - 4A±fc(0) = 4/3^(1) . 

Here, the so-called "generalized momenta" P ik - do not depend on velocities, which appear 
only on the right hand side of the equation. Explicit calculation leads to the result 
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= Aa[AT i±k + \BT- kA2 + \C m - k T^ ±m + \{B + C) ni T ±± - k ] . 

This system of equations can be decomposed into irreducible parts with respect to the 
group of three-dimensional rotations in Eo: 

P±k = *_ui/J - 2aCT m mk = 2a(2A + B + C)T ±±k , 
Pi = *ik/ J ~ 2aBT m = 2a(2A - B)T{^ , 
P^^/J = M^ + B)T{ ±k , 
P m m = ir m rn/J = 2a{2A + B + 3C)T rT 



-rm 

Lm j 



where X|| = Xhm, XT^ = Xnj^ — nj k X n n /3. Taking now into account the special choice of 
parameters adopted in equation ( |2.2| ), we recognize here two sets of relations: the first set 
represents extra primary constraints, 

P±k = 7T±fe/ J + 2aT n l mk « , 

Pi=^i/J-aT 1 ^^0, (3.4a) 

usually called if- constraints , while the second set gives nonsingular equations, 

P T=^/J = 2aT{ ± - k , 
P^ m = ^ m /J = -4aT^ ±m , (3.4b) 

which can be solved for velocities. 

Further calculations are greatly simplified by observing that both sets of extra constraints 
( p.4a| ) can be represented in a unified manner as 

<fe = n& - ir M + aV a B% , B% = effi^ . (3.5) 

This is seen from the fact that relations (|3.4a|) can be equivalently written as 

^ik ~ n ki ~ 2a J (TjJk — n iT m rhk + n kT m ml) = 2aV aH® k , 



where the last equality follows from (|A1|) , and the identity 2H^ k = —B^ k . 

2. Having found all the primary constraints, we now proceed to find the canonical 
Hamiltonian density ||22|| : 

H c = TT i a b i a + \>K i fA i i tt -bC. 
The velocities b l a and A^ a can be calculated from the relations defining T l 0a and R^oa'- 

T 0a = d b\ + A l m0 b m a - d a b\ - A l ma b m = NT\ a + N^T Pa , 
R ij 0a = d A l \ + A l m0 A m i a - d a A ij a - A l ma A m \ = NR ij ± a + . 

After a simple algebra we find that the canonical Hamiltonian can be written as a linear 
function of unphysical variables (b k o, A^ ), up to a 3-divergence, 

H c = NH± + N a H a - \A i] H i3 + d a D a , (3.6a) 
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where 

Hij = 2'K[ i l3 b j ]p + V a ^ij a , 

'Ha = ^k^T k aP - b k a Vp^k^ + fa/ R lj a/3 , 

H± = {t k T\ k - JC T - n k V^tt/) - JX^RV^ , 

D a =b k 7r k a + lA i i ir ij a . (3.6b) 

Here, TCij and 7i a are purely kinematical terms whose form does not depend on the choice 
of the Lagrangian, and T~C± is the only dynamical part. 

The explicit form of 7i± can be obtained by eliminating "velocities" T i±k with the help 
of the relations defining momenta ir ik . To do that we first rewrite the first two terms of 7i± 
in the form 

K lk T i± - k - JC T = \ JP %kr ^n_k - J£t{T) , 
where T ik i = T ik j. Then, taking into account the constraints (|3.4a|) one finds that T ±±k 
and TAt are absent from TC±, whereupon the relations ( |3.4b| ) can be used to eliminate the 
remaining "velocities" Tj^t and T m ± m , leading directly to 

n± = (|P| - JC T (f) - n fc V^7r/) - JXif^rnn , (3.7a) 

where 

P 2 1 / (Ik) 1 m Jh 

T 2aJ V ( } 2 

7> (rp\ „ / Irp _n-\rnnk < It _rpnmk rpfh _rp nk\ /q *7V^ 

'-'TK- 1 ) — « yz 1 rank 1 +2 1 rrink 1 ~ 1 fhk 1 n J ■ {6. ID) 

The general Hamiltonian dynamics of the system is described by the total Hamiltonian, 
which is given as 

+ §u < W', (3.8) 

where u's are, at this stage, arbitrary Hamiltonian multipliers. 

Although the torsion components T ±±k and T-^r are absent from the canonical Hamilto- 
nian, they re-appear in the total Hamiltonian as the non-dynamical Hamiltonian multipliers. 
Indeed, the Hamiltonian field equations for b k a imply (Appendix B) 

NT ±±k = u ±k , NT\ ±k = vt k . (3.9) 

The presence of non-dynamical torsion components does not imply that GRii is an inconsis- 
tent theory [Q, as it has very clear interpretation via the gauge structure of the theory: it 
is related to the existence of additional first-class constraints fau, as we shall see in Sec. V. 



IV. CONSISTENCY CONDITIONS 

The consistency of the theory requires that the constraints do not change during the 
time evolution of the system governed by the total Hamiltonian: 

X {x) = j t (f>{x) = |rfV{0,K)^O, 

where {A, B'} denotes the Poisson bracket (PB) of two variables A(x) and B(x'), and 
x° = (x')°. The integration sign will be often omitted for simplicity. 
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A. Consistency conditions of primary constraints 

Having found the form of primary constraints in GRy, displayed in Eqs. (RTTT), ( p.2| ), 
73]) and ( |3.5| ), we now consider the requirements for their consistency. 
Since the canonical Hamiltonian is linear in unphysical variables (b l o, A^q), the consis- 
tency conditions of the sure primary constraints are given by 

x± = n±^o, X a = n a ^o, X y = ^«o. (4.1) 

By noting that the components ir^ a p in Eq.( |3.2|) have vanishing PBs with all primary 
constraints, we easily obtain 

X ij af3=R ij af3~0- (4.2a) 

On the other hand, the consistency of 7r* J oa implies 

tf j ofi = u*f i -N a ff* a f l K i O => u ij p*0. (4.2b) 

The dynamical meaning of the last condition can be seen more clearly if we note that the 
equation of motion for A % i p, 8qA 1 ^ p = {A lJ p, H c } + u^p, can be transformed into the form 

RV p&u ij p. (4.3) 

Hence, Eqs. ( |4.2j ) tell us that all the components of the curvature tensor weakly vanish, as 
one could have expected. 
Using the PB relation 

<M = a( Vik B% + Vjk B° a )S - (*0 , 
the consistency condition for 0^° takes the form 

Xl3 a = {<Pij a , U c } + a( Ul s B° s f + Uj 'B°?) - Abu^ « . 
It can be used to determine Uij 0a : 

Uij 0a = ^ S B° S « + Uj s B%) + , fiy«* = 1{ V, H c } , (4.4) 

where {iiij a ,T-C c } is calculated in Appendix C. The first part of Uij 0a contains u k i and gives 
an additional contribution to u kl <p k i, so that the replacement of this result into Ht leads 
effectively to 

4> kl = <f> kl - ^(7r k s oaB° s ? + TTfoaB^) . (4.5) 

Note that {fa^, <f)' kl } = 0. 

The most complicated consistency conditions are those for the tetrad constraints <pij (or, 
equivalently, <f>ij). First we note that their PB algebra has the form 

{0tj, <P'mn} = (Vim&nj + Vjm&n)^ ~ (mn) , (4.6) 

and that the term \u l i a (f)ij a in TCt can be discarded according to ( |4.2b|) . Then, after showing 
that the Poisson bracket {<fiij,T~l r c } vanishes weakly, we will be able to conclude that the 
consistency condition for 0^ is automatically fulfilled: 

Xa : = {^-,^}«{^,K}«0. (4.7) 
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B. The consistency condition of <pij 

In order to simplify the derivation of the consistency condition for (j)^ we rewrite this 
constraint in the form 

4*ij Fij , 

Fij = V a (V - aB°?) = VJi^ . (4.8) 

General arguments in PGT, related to the local Lorentz symmetry of the theory, imply 
that the constraint TCij is of the first class p2 |. This is also clear from the PB algebra of 



constraints, discussed in the next section. As a consequence, the consistency of 0^ follows 
from the consistency of F^. 

We are now going to show that {Fij,H' T } ~ 0. First, we note that 

{%0' fc J = O. (4.9) 

Then, using the results (CAa) we obtain 

{Fij, H' kl } = (rjikFij + VjkFu) 5 - (kl) , 

{F^H'p} = -VpfaiS) + [V a , Vp](U l3 a S) . (4.10a) 
while (CAb) leads to 

{Fj, H'J = \V p [(n4 jk - n j( p ik )h^5} + ±[V a , Vp}(M^5) , 

NMf = 2aHf - 46A,/ Q + N a (U t / - 0,/) - A^IL/ " <f>ij a ) , (4.10b) 
Hence, the consistency condition of Fy, and consequently of (f>ij, is automatically satisfied. 

C. Consistency conditions of secondary constraints 

In the process of investigating the consistency of primary constraints in GRy, we obtained 
secondary constraints ( [i~TD and (|4.2a|) . 



Consider, first, the consistency condition of the secondary constraint R 13 a p- Since R lJ a 
depends only on A tJ a , one can express dR 13 ' a p/dt in terms of dA ZJ a /dt, use the equation of 
motion ( |4.3| ) for A % 3 a , and rewrite the result in the form 

V R i3 aP w V a u i3 p - Vpu i3 a . 

Hence, the consistency condition for R lJ a p is identically satisfied. 

The above relation has a very interesting geometric interpretation. Indeed, using equa- 
tion ( |4.3| ) we see that it is a weak consequence of the second Bianchi identity. 

General arguments in PGT show that the secondary constraints r H.ij, r H, a an d 7~t± are 
related to Poincare gauge symmetry |p2| , p4"|| . Consequently, they are of the first class, and 
their consistency conditions are automatically satisfied. This will be explicitly seen in the 
next section, from the form of their PB algebra. 

Finally, at the end of the consistency procedure, we give the final expression for the total 
Hamiltonian: 
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H T =H' + TiVi° + + \Uij a ^ ij aP + \^fa , 

H'= H e + 1^^0/3. (4.11) 

The multipliers u 1 q, u l \, Uij a/3 and u lk remained arbitrary functions of time, hence we expect 
that the related constraints 7Tj°, ir^ , a/3 and fa are of the first class. 

V. THE ALGEBRA OF CONSTRAINTS 

In the previous analysis we found that GRn is characterized by the following set of 
constraints: 

primary: 7T;°, ir^ , fa, it 13 a/3 , fa a , n l: > p ; 
secondary: 7ij_ , H a , Hij , R v a p . 

It is simple to see that fa 01 and n^op are second class constraints. They can and will be 
used as strong equalities to eliminate Xij 0a and tr^op from the theory and simplify further 
exposition. In particular, the term with the determined multiplier u^ a in the total Hamil- 
tonian can be now neglected, and fa reduces to fa, equation (|4.5|) . Due to a simple form 
of these second class constraints, the related Dirac brackets have the form of PBs in the 
phase space of the remaining variables. All the remaining constraints are of the first class, 
as follows from their PB algebra. 



Since the kinematical constraints TCij,TC a ,TC± have the same general form as in p2| . |24 
their algebra remains the same: 

{Hij,H' kl } = {r)i k Hij + r] jk Hu)5 - (kl) , 
{Hij,H' a } = , 

{H a , H'p] = (H' a dp + Hpd Q - \R\pH i3 )b . (5.1a) 

As a consequence of R 1 ^ a p ~ 0, the last term in {Ti^T-Cp} is quadratic in constraints. 
The brackets involving Ti,± are found to have the form 

{Hij,TC' ± } = 0, 
{H a , H' ± } = H±d a 5 , 

{Kl, H' ± } = -C9 aP n a + 3 g' aP H' a )dp6 . (5.1b) 

The first two brackets are most easily verified by taking into account that 7i± can be 
written in the form Ti,± = J f(£, A ) — n k, V a iik a , where / is a Lorentz scalar formed from 



variables £ = (T^n, R lJ m,n,^i / J,^ij / J,n ), as shown in fl24[. The second bracket is 



obtained from the general formula based on the chain rule for PBs, 



a 



which explains why the second term is absent in (|5.1b|) . 

The last and most important bracket {Ti^,Ti' ± } is evaluated using the chain rule and 
keeping only those terms that contain d a 5 (terms proportional to 5 do not have the correct 
symmetry under x <-> x', hence they cancel each other). 
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In the next step we want to extend the above algebra by adding </>y = <pij. The relevant 
PBs involving faj are given by 

{4>ij, 4>m} = {vik<Pij + Vjk4>u) - (kl) , 

{0tj, Wju} = (Vik4>ij + Vjk4>u)$ - (kl) , 
fa,, wy = - [V a , VfAQlv-S) , 

= -\V p [{n i( j> jk - n j( j> ik )h^8] - §[V a , Vp](M^8) . (5.2) 

Finally, we display the non-vanishing PBs involving R lJ ' a p and 7r l3 a p: 

{R ij a p,H' kl } = (W*? + 6{R\ a/3 )8 - (kl) , 
{R i3 a p,H^} = V a (R i3 lP 5) - (a/3) , 

K J a/3 ,^l} = 4J^ a/3 . (5.3) 

Thus, all constraints except 0jj a and if" 3 op are of the first class. The fact that <f>ij is 
first class is of particular importance for the consistent interpretation of the non-dynamical 
torsion components, as noted at the end of Sec. III. 



VI. EXTRA GAUGE SYMMETRIES 

The presence of arbitrary multipliers in the total Hamiltonian is related to the existence 
of gauge symmetries in the theory. The general method of constructing the generators of such 
symmetries has been given by Castellani p5| . If we limit ourselves to gauge transformations 
given in terms of arbitrary parameters e(t) and their first time derivative s(t), which is 
sufficient for the present analysis, the gauge generators take the form 

G = J d 3 x[e(t)G {0) + e(t)G {1) ] , (6. la) 

where and G^ are phase space functions determined by the conditions 

G = Cppc i 
G® + {GV),Ht}=Cp F c, 

{G ( -°\H t }=C pfc , (6.1b) 

and Cppc denotes primary first class (PFC) constraint. 

The Poincare gauge symmetry is present in our formulation of GRn by construction, and 
the related gauge generator is based on the sure constraints 71^,71^° and H±,TC a ,TCij p6 



Here, we shall focus our attention on extra gauge symmetries based on 7T l3 a p, and R l3 a p- 

A. Extra gauge symmetry, 1 

Starting with ix %3 a p as G^ in (|6.1b|) we find that the related gauge generator is given by 



G = d 3 x 



f (Vo% a/ V Q /3 + H- a/3 (-4W^ + (6/6)7^)1 . (6.2) 
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The only nontrivial gauge transformations 5qX = {X, G} are 

5 (b\ t ^) = V (fcii a/3 ) , 

5 n tJ a = AV p {be l3 ^) . (6.3) 

To see the meaning of these transformations, consider the Hamiltonian equation for the 
variable Hij a = it^ 01 — aBf". Introducing K°f = AbX^ 13 — aB°f and using the results of 
Appendix C we obtain the equation 

VoIV-V^f = 0, (6.4) 

which is the Hamiltonian analogue of ( |2.4b| ). The application of the above gauge transfor- 
mation to this equation yields 

(V V^-V^V )(4fe^) = 0. 

The invariance follows from the fact that the left hand side vanishes in Weitzenbock space, 
where R^op = 0. 



Starting with G^' 



where (Appendix D) 



B. Extra gauge symmetry, 2 

bij in ( 6.1b ), one finds that the gauge generator has the form 



Gij = J cPxQ&gW + \e ij Gf)\ 



G 



(o) 



1 1 



r> s TS&ft i / A n „ s i as „n \ T^afB s TS&fi 



(6.5a) 



(ij) . (6.5b) 



The corresponding gauge transformations are: 
5 Q b k a = £ k s b s a , M° a = 0, 

Ab5 X^ = \Ei n k^ + ATo(e m n K$ + e^K*)] - (ij) , 
8 mj a = WfB^) + V^Cf )] - {ij) , 



(6.6) 



and similarly for other variables. 



Consider, again, equation (|6.4| ). Using 5oIlij a = V p(ei n K^) — (ij), we easily obtain 



^(v n !J Q ) = v ftn ! 



nj J 



(ij) 



VpWKtf + e?K% + Ai' (e»K$ + e?K%)\ - (ij) 



where we made use of 



0/3 



0. On the other hand, 



SoKg = [ei n K? + efkg + A l %(e s n K^ + efIC£)\ - (ij) 



and we see that equation ( |6.4|) is gauge invariant. 
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VII. CONCLUDING REMARKS 



The investigation of the Hamiltonian structure of the teleparallel formulation of GR 



presented here is based on Dirac's general method for constrained dynamical systems |20 |. 

To complete our results, we now discuss how the physical degrees of freedom of GRy 
are counted. After the elimination of Xij 0a and ir^o a , the reduced phase space is spanned 
by the 40 + 18 field components (b\, M , X^ 13 ), and the same number of momenta. The 
primary first class constraints 7r lJ ' a p diminish the number of independent variables for 2 x 18, 
leaving us with the phase space containing effectively 2 x 40 components. Before going 
on, we wish to clarify the counting of constraints R 1 ^ a p ~ 0. Note that here we have 
formally 18 equations, but they represent only 12 independent conditions on A % \. Indeed, 
starting with the simplest solution A^ a = of R^ a p(A) = 0, one can construct a new, 
Lorentz rotated solution A^ a (A) = A l kd a A^ k , containing 6 arbitrary parameters A lk ||, so 
that the number of independent conditions on A^ a is 18 — 6 = 12. Continuing now the 
counting, we find 20 sure first class constraints [ten primary (fti , 7Ty°), and ten secondary 
(H±, H a , Wjj)], and 6 + 12 = 18 additional first class constraints 0^ and -R* J a /3, which leaves 
us with 2x40 — 2x 38 = 4 physical degrees of freedom, corresponding to massless graviton. 

We found two types of extra gauge symmetries in the PGT formulation of GRy. The 
first type is related to the primary constraints 7r u ai g. The related gauge transformations do 
not act on b l ^, hence they are irrelevant for the structure of the first field equation ( |2.4a ). 
On the other hand, the gauge symmetry acts nontrivially on Lagrange multipliers. If we 
recall that the only role of the second field equation ( |2.4b| ) is to determine these multipliers 
||, it becomes clear that this cannot be done uniquely without fixing the gauge. 

The second type of extra gauge symmetry originates from the tetrad constraints <fiy. 
We note that Nester |L4| derived these constraints in the form ([3.4a ), in his analysis of the 
positivity of energy in the teleparallel form of GRm. Their existence may be interpreted as a 
consequence of the fact that the velocities contained in T ±1 j, and T^r appear at most linear in 
the Lagrangian [13] and, consequently, remain arbitrary functions of time. The phenomenon 
that some velocities are dynamically undetermined is quite usual for constrained dynamical 
systems PU[. For the related initial value problem to be well defined, these undetermined 



velocities should be removed from the set of dynamical velocities [T^ . 



The role of this symmetry is very clearly seen if we observe that the teleparallel geometry 
can be also formulated as the translational gauge theory, where local Lorentz symmetry is in 
general absent |5]||. However, for the special choice of parameters corresponding to GRy one 
finds that cf)^ is an additional first class constraint, which generates local Lorentz symmetry 
as an extra gauge symmetry [Q. This also clarifies the form ( |3.5|) of <pij, which is seen to 



"imitate" Hij in the tetrad sector. 

Maluf [21j studied GRy by imposing the time gauge at the Lagrangian level. His argu- 
ments concerning the necessity of the time gauge in the canonical formalism are conceptually 
misleading: this gauge (as well as any other gauge) may be useful, but certainly not essen- 
tial PUj. After fixing the time gauge, he found the Hamiltonian and derived the constraint 
corresponding to our 0^. [Eq. (25) in his paper], while is missed. Moreover, Maluf was 
not able to calculate the constraint algebra unless imposing another gauge condition. His 
constraint algebra [Eqs. (30)-(34)] does not agree with our results, which might be a conse- 
quence of the adopted gauge conditions. All this makes this analysis of the gauge structure 
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of GRy rather unclear. 

The results obtained in this paper refer to non-interacting GR|| , and can be used to define 
and analyze the gravitational energy and other conserved quantities [0,0 • Interaction with 
matter fields may be included in a straightforward manner [0,^1 . Studying consistency 
requirements imposed by extra gauge symmetries on the matter sector will tell us more 
about the existence and nature of consistent couplings ]n| . 
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APPENDIX A: SOME GEOMETRIC IDENTITIES IN T 4 



We begin with a simple but technically important identity 



hi/ 



bh k »(T k tl - b k T 3 + 5jTi) = -4b/3 {ij f/a . 



H% = bihfhj 1 



h/hi 



(Al) 



which implies V ^{b(5 {ij f ) = for R\ v = 0. 

In Riemann-Cartan space U4 the Lorentz connection can be expressed in the form A = 
A+K, where A is Levi-Civita connection and K the contortion. Substituting this expression 
into the definition of the curvature tensor R l ^ u (A), we obtain the basic identity: 



(A2) 



where V = V(A) is Riemannian covariant derivative. Then, multiplying this relation with 
Hg/2 and using V' M #|f = 0, we find 



bR(A) = bR(A) + b(\T tjk T^ k + \T l]k T> lk - T k T k ) + 2d,(bK») , 

where K» = K^ n n = — T M . 

Now, if we write equation (A2) in an equivalent form, 

#V(A) = #V( A ) + [V,K 1 \ - K\^K S \ - (ja 

and multiply it with H k J /2, we obtain the result 

bR\{A) = bR\{A) + 
which can be written as 



(A3) 



n jk 1 



abR ik (A) = abR tk (A) + 2V M (fy3 v/c ) + 2bf3 mn k T 



yik 



krpmni 



-b(3 imn T k mn - v lk ad^bT») - 4V M (6/5 



ik]fi\ 



The last term on the right hand side vanishes for R^ pv(A) = 0. In that case we find 
2ab[R ik (A) - y h R(A)] = -4V M (6/3 vfc ) - Ab(3 mn k T mni + 2bp imn T k mn + if k bL T . 



(A4) 



(A5) 
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APPENDIX B: UNPHYSICAL TORSION COMPONENTS 



In this Appendix we show that the unphysical torsion components T ±±k and 7V}j can be 
expressed in terms of the Hamiltonian multipliers u ki . 
Using the PB relations 

{b\,H' kl } = SlkJ - (kl) , 

{b\, H'p} = (y a b i p )6 - b l p & a 5 = T p J + V a (6» , 

{b\,H' ± } = ± (b ma P^ - ^ai^m) 5 + V Q (n l 5) , (Bl) 

one easily finds that the Hamiltonian equation for b k a can be written in the form 
V b\ = Va&'o + N^T Pa 

+Nha ^J ^ ~ + bk « + uiS ) • ( B2 ) 

As a consequence, 

u ±k = NT ±±k ^ = jyyUfc _ (33) 



APPENDIX C: CONSISTENCY CONDITIONS 



We collect here several technical relations which simplify the derivation of the consistency 



conditions for 0^° and F^. 



1. The term {irij a ,H c } in the consistency condition for the primary constraint <pij a is 
calculated using the relations 



{-Kij a ,H' kl } = (rjiknij 01 + r) jk 7r u a )5 - (kl) , 



(Cla) 



and 



K<\ H' ± } = -4V /3 ( JX t / a 5) + AV p [J(N a \W - N^^Sl , 
- 8J(3 m (0) + 2aJn {i T ±]]k ] h la 5 , 
+ {niiif k - nj-K ik )h ka 5 . 

Using the identity 8J(3 mk {0)h ha = ae^ n T m lfi n n , we obtain 

{^ij a , H c ] = -(Ai S 7T sj a + Adonis 01 ) - iV a (7T- - TTjl) 

-4V^(6A,/ Q ) - aNe^ n T m lP n n + A^7r ( ^ - n^h** 



(Clb) 



(C2) 
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2. In order to calculate the Poisson brackets between and the Hamiltonian con- 
straints, we also need the following relations: 

{B^, H' kl } = ( Vik B?° + r] jk B^)S - (kl) , 

{B™, H'p] = Vp(B%8) + tftfS&J , (C3a) 

and 

W, H' ± } = Un iHTk) - n^ m )h ka 5 + 2eJ^V 7 (n n <5) . (C3b) 
Combining (Cla) and (C3a) we find 

{IV, H' kl } = (Vik~ttij a + V^") 6 - (kl) , 

{IV, %} = -5^5 - <^V 7 (IV<5) + V^(n ii a 5) , (C4a) 

which implies ( [4.10a ). 

Similarly, combining (Clb) and (C3b), and using the identity 

Bf = 2Nef2X nn - ( NaB if - N " B "D > 

we obtain 

{IV, tt^} = f (n^ - n^hH + V p {M«f5) , 

NMf = 2aH?f + 4&A^ + iV Q (IL/ - 0,/) - A^(IV - V) , (C4b) 
which implies ( |4.10b|) . 



APPENDIX D: EXTRA GAUGE GENERATORS 



In this Appendix we derive the form of the gauge generator ( Of) . We start with G 



(i) 
■i 



'13 



in ( |6.1b| ). In order to find the form of the accompanying component G\j\ we use the PB 
algebra given in Eq. (|5.2| ), and calculate 



{<fa,H T } = -l\y a ,Vp]K 



a/3 
ij > 



K°f = 2aHff + 46A 



a/3 
ij , 



where terms proportional to <fiij a are discarded. The second condition in ( |6.1b| ) implies 

(Dla) 



where #y is a primary FC constraint. The third condition in (|6.1b|) can be written in the 
form 



\[Ri S apK^j + Ri S a/3 Kgf — (ij)] + 9ij = Cpfc , 



where X = {X, H?}- Then, using the relations Vq-R. 



i a(3 



and 7T l 



ij - 2 46 



(Ai n 7l n S al 3 + A S n0 7li n a) 3) Kgj — -KiafiK, 



' a p = 4bR ls a p we obtain 



(ij) . (Dlb) 
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